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We study theoretically the bulk modulus (inverse of the compressibility) of a suspension of charged 
objects (macro-ions), making use of a cell model to account for the finite density of macro-ions. The 
O ■ diffuse layer of charged micro-species around a macro-ion is described by a generic local density 

D ' functional theory. Within this general framework, we obtain the condition for a positive bulk 

' modulus, which is fulfilled by several proposals made in the literature and rules out the possibility 

of a critical point. We show that a sufficient condition for a positive compressibility also ensures 
that the same theory produces repulsive effective pair potentials. 



I. INTRODUCTION 



Macromolecules soluble in aqueous suspensions usually acquire an electric charge. Such systems are of considerable 
theoretical and experimental interest. Examples include proteins in living cells, dispersion paints or superabsorbants. 
■ Their theoretical description is however a tremendous task, and operational approximate treatments are needed. 
\ In this work, we consider a cell model to account for the macro-ions correlations^, in conjunction with a local 
P5 ' density functional theory to describe the inhomogeneous electric double-layer around a macro- ion^''^. This framework 
O " encompasses the standard mean-field Poisson-Boltzmann (PB) theory, but also more recent approaches proposed to 
, account for micro- ions excluded volume^, more general non electrostatic effects^ or to incorporate correlations between 

the screening micro-ions®'^, that are neglected within PB theory. The cell model description may be considered as 
one of the simplest starting points and provides an important benchmark against which experiments and more refined 
^ theories are tested. The resulting differential equations are however highly non-linear, and even with the simplest of 
\l I the approaches under consideration here (PB), can only be solved analytically in 1 or 2 dimensions without added 
salt (i.e. for the situation of a flat double-layer in a confining slab, or for that of a rod-like macro-ion of infinite 
length enclosed in a concentric confining cylinder^). The linearized version of the above problem has therefore always 
. been an important alternative, but is not free of internal inconsistencies. In particular, the linearized PB theory may 
' lead to negative osmotic pressures^' '^'^ and negative bulk modulus (inverse of the compressibility)^*^, whereas within 
, the original non-linear PB theory, it is straightforward to show that the osmotic pressure is necessarily positive. For 
• the bulk modulus, the situation is less clear: although this quantity is found positive in the numerical solution of PB 
equation, it seems that no formal proof exists concerning its sign. In this article, we derive such a proof and obtain 
the conditions under which a general local density functional theory leads to a positive compressibility within the cell 
, model. 

. The article is organized as follows. The density functional theory formalism is presented in section II where a few 
Q ' useful identities are derived. The bulk modulus is then computed in section III and cast in a form where a sufficient 
^ . condition for its positivity clearly appears. Conclusions are drawned in the the final section. 



II. GENERAL FORMALISM 
A. Density functional theory 



We consider a z/-dimensional colloid with spherical symmetry confined in its (concentric) Wigner-Seitz cell [v — 2 
for a cylindrical colloid, v ^ 2) for a colloidal sphere). The cell is a spherical region TZ of radius R and volume 
V = S,jR'^ /v {Si, is the area of the unit radius z^-dimensional sphere, ^2 = 2tt, S3 = 47r). The colloid (of fixed uniform 
surface charge density) is immersed in an electrolyte solution with several different species of ions with charges {(/a}, 
and the local density of the species a is denoted na{r). The system inside the cell is globally neutral. In what follows, 
we explicitly consider the semi grand-canonical situation where the macro-ions suspension is in osmotic equilibrium 
with a salt reservoir through a semi-permeable membrane (permeable to micro-species only). We therefore consider 
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the grand potential O {n^}, which is a functional of micro- ion densities that we write 

n{nc.}^ [ a;({na(r)})dr+^ / p(r)V(r) dr + A / p(r) dr, (2.1) 



n 2 Jn j-fi 



where 

'°('') = X] 9a^^a(^') + Pcoi(r) (2.2) 

is the total charge density including that of the colloid (denoted Pcoi); A is a Lagrange multiplier to ensure global 
electroneutrality {dfl/dX = J p = 0) and ip{r) is the electrostatic potential. Within PB theory, the local part of the 
fimctional (2.1) embodied in uj is entirely of entropic origin. In other words, electrostatic interactions are taken into 
account at mean-field level only, through the term J p'il'/2. This feature is shared by the theories of Refs. [4,5], but 
is not the case of the formalism put forward by Barbosa et al.^'"^ , where the term uu contains electrostatic corrections 
to the mean-field Coulomb contribution. By comparison with the results of Molecular Dynamics simulations, this 
framework was shown to capture important correlations missed by the mean- field Poisson-Boltzmann^. 

Within the standard Poisson-Boltzmann theory, the micro-ions are treated as an ideal gas of charged parti- 
cles reacting to the mean electrostatic potential. The free energy density therefore reads /({«■«}) = /id ({«■«}) = 
(3^^ na{ln[A^na] — 1} where A is an irrelevant length scale and /3 = l/{kBT) the inverse temperature. The grand 
potential functional describing the osmotic equilibrium with a salt reservoir is thus = J f — l^a J ''^a ^ J pV-'/2 + 
Xjp, where the chemical potentials p.'^ = 9/id({n^})/9n^ = In(A^n^) are defined from the bulk densities of 
micro-species in the reservoir. The local grand potential density finally takes the form 

a;({n„(r)}) = /({n«(r)}) - ^M^n„(r) = ^n„(r) (in^ - l) . (2.3) 

As mentioned in the introduction, other theories may be described by our formalism, such as those proposed to account 
for stcric effects'*, more general non-electrostatic interactions^ or to go beyond mean field and include correlations®'^. 

The mean- field Coulomb term ^ J p(r)f/;(r) dr in Eq. (2.1) is actually quadratic in the densities if one writes the 
potential in term of the set {n^}. Introducing the proper Green's function G(r, r') for the electrostatic problem in 
the region IZ, it is always possible to write the electrostatic potential in the form (see appendix A) 

V'(r) = / p(r')G(r,r')dr', (2.4) 



JTZ 

with the cell boundary chosen as the potential refcrcnc(\ ?/'(/?,) = 0. 

The variational problem defined by the functional (2.1) leads to the stationary condition 

'"y^^^ ..(^(r) + A) = 0, (2.5) 

which, in the Poisson-Boltzmann theory, reduces to 

n„(r) =<e-^«°['^('-)+^l. (2.6) 

Since the potential vanishes at the boundary of the cell, it can be seen in this equation that the Lagrange multiplier 
A coincides with the so-called Donnan potential (potential drop across the interface of the cell and the bulk reservoir, 
which is formally at a potential —A). 

Finally, discriminating between the optimal density profiles fulfilling Eq. (2.5) and the generic arguments of the 
grand potential functional leads to unnecessary heavy notations and will not be useful for the subsequent analysis. 



B. A few useful identities 



Consider a functional Q of the density profiles na{r) having an explicit dependence on the volume V of the cell 
and the potential drop A (later, we shall be interested in Q = O). The total derivative dQ/dV of Q when the volume 
of the cell varies is^° 

dQ _dQ ^ /■ V- duaiv) dQdX 



dV dV Jn'^^nair) dV dX dV 
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where the first term is tlie variation due to the explicit dependence of Q on V, the second is due to the variation of the 
density profiles when the volume changes, and the last is due to the variation of the potential drop with the volume. 
Computing dna{'r)/dV and dX/dV would require to solve the variational problem (2.5) subject to the neutrality 
condition for a cell of volume V + dV and for a cell of volume V, before computing the difference of the two solutions. 
However this will not be necessary for our purposes. 

We will be interested most of the time in quantities defined by a local density. In general the partial derivative with 
respect to the volume (explicit dependence) of such a quantity may be computed by means of the dilatation method: 



d_ 
dV 



g{v) dr 



n 



d_ 
dV 



r-^dr 
or 



V I g{V''/''r)dr 

. Jiz 

= y f 5(r)dr+-i- / g{r)vdS-^ [ ug{v)dv 



)r • dS. 



(2.8) 
(2.9) 



We first made a change of variable r = V^/'^r to show explicitly the volume dependence of the integral. After computing 
the derivatives and returning to the unsealed variable r we made an integration by parts {v is the dimensionality of 
the cell = 2, 3). For the spherical isotropic geometry we are interested in, this result reduces to 



d_ 

dV 



g(r) dr 



n 



9{R)- 



(2.10) 



We will need also to compute the derivative of some terms given by a double integral. Following the same steps as 
above, this derivative is given by 



d_ 
W 



[ g{r,r')drdr' = [ g{R,r')dr'+ [ g{r,R)dr. 
Jn J Jtz Jn 



(2.11) 



The key assumption here is that both quantities / g{r, R) dr and / ^(R, r) dr only depend on the modulus R= |R|. 

Another important point to notice is that the volume total derivative is taken at constant total electric charge because 
the system is neutral. Applying equation (2.7) to qoi^a = / Ylia la^ai'^) '^^ then, gives 



dN^ 



= P(^)+XE-^^r. 



(2.12) 



III. THE COMPRESSIBILITY 



A. Derivation from the grand potential 

Applying formula (2.7) to the grand potential once gives minus the pressure. However, as noticed in Ref. [10], the 
second and third terms vanish at the solution of the variational problem. Making use of equations (2.10) and (2.11), 

we consequently have 



— = ^{MR)}) + p{RmR) + >] 

= w{{noc{R)}) + Xp{R). 



(3.1a) 
(3.1b) 



The last equality follows from '>p{R) = [see the property (AlO)] and is valid in general for arbitrary isotropic 
densities even if they differ from the optimal ones solving the variational problem. This will be used later. Note that 
the symmetry property of the Green's function [G(r, r') = G(r', r)] is an important ingredient in obtaining the above 
equations^^. 

Using the stationary condition (2.5) one finds the usual result^'^'^° 



P = ^na(-R) 



dna{R) 



-a;(K(i?)}) 



(3.2) 
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which reduces to p = A;sT^^na(-R) within PB theory (see Ref. [13] for a general derivation of this result, valid 
beyond PB). 

When computing the second total derivative with respect to the volume of O to obtain the bulk modulus, one 
should not disregard the second and third terms of equation (2.7) to early. 



dV2 
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+2 
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drig (r) 
dV 



dr 



dVdX dV 



„ f ^ rsri] dna(r) dX 

"1? 



dX 



n 2 



dV dV 



dV 



(3.3a) 
(3.3b) 

(3.3c) 
(3.3d) 

(3.3e) 

(3.3f) 



At the solution of the variational problem most of these terms vanish, as it will be shown below. The first term (3.3a) 
is obtained taking the partial derivative of (3.1b) with respect to V 



d_ 
W 



dna{R) 
dR 



duj 



dna{R) 



= 0. 



(3.4) 



We have used the variational equation (2.5) at r = R e dTZ and the fact that i^{R) = 0. The second term (3.3b) is 
obtained replacing the expression (3.1b) for dfl/dV into (3.3b) 



dna{R) 



+ qaX 



dna(R) 
dV 







(3.5) 



where wc have once more used equation (2.5) at r = R S dTZ. 
The third term (3.3c) is equal to 



and with the fifth term (3.3e) 



d 



n 



^+<Z.(^(r) + A) 



dna{r) dX 
dV 



dn„(r) dX , ^ dX 
,^-c/r = -2p(i?)- 



(3.6) 



(3.7) 



gives a vanishing contribution. Wc have used in the preceding equation the relation (2.12). The last term (3.3f) is 
zero since in linear on A. Finally the inverse compressibility may be cast in the form 



From equation (2.1) we have 



dp 
dV 



yP.=^[ ErT^^%^^^r.r'. (3.8) 

"•^ i-R2 onc((r)c)n-y(r') dV dV 



The second term, the Coulomb contribution, when replaced into equation (3.8) gives 

, dna{r) dn j{r') 



dV dV 



■ drdr' = — ^ [ (/>(r)A0(r) dr 



-R^-^mdnm 



1 



|V(^(r)|2 dr 



n 



(3.9) 



(3.10a) 



(3.10b) 
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where we have defined the "electric potential" created by the charge variation 

0(r)= / G(r,r')^fc- 



dV 



(3.11) 



and performed an integration by parts. The boundary term in (3.10b) vanishes because, from equation (AlO), one 
has 



m = / E 



g„^^G(R,r')dr' = 0, 



thus showing that the Coulomb contribution term is always positive and 

d'^uj dna{r) dn^(r) V 



X 



n~ dnadn-y dV dV 



+ ^ / Vr / G(r,r')E9< 



' dV 



dr. 



(3.12) 



(3.13) 



Discussion 



From Eq. (3.13), the positive definiteness of the integral operator whose kernel is defined by 



dnadn-y 



6{v - r') 



(3.14) 



ensures that the compressibility is positive (this is a sufficient but not necessary condition). This is the case for 
Poisson-Boltzmann theory where u) is simply the ideal gas grand potential density, and it may be checked that it also 
holds for the theories presented in Refs. [4-6]. 

More generally, in any well constructed approximate theory for the colloid based on density functional for the 

grand potential of the form (2.1) (that is a local density term plus an "interaction term" given by the mean- field 
Coulomb electrostatic energy), the solution of the variational problem should be a minimum, i. e. the quadratic form 



5na{v)5n~j{Y') 



(3.15) 



should be positive definite to ensure thermodynamic stability. From expression (3.8), it then follows that the com- 
pressibility will be always positive in such a theory. 

Equation (3.8) for the compressibility can be seen as a generalization for non-uniform fiuids of the compressibility 
sum rule for uniform fluids^'* 



/3("X) 



c(2)(r)dr 



(3.16) 



written here for a one-component system. In this relation, S-'^\r) is the direct correlation function defined, in the 
more general situation of a mixture, from 



(5na(r)5n^(r') 5na{r)6n-y{v') 



-kBTc'^^]{v,v') 



(3.17) 



where J^i^ is the ideal gas contribution to the free energy functional. When one replaces (3.17) into (3.8) for a uniform 
fluid [ca-^(r, r') = Ca-ydr — r'|)], and knowing that for a uniform fluid dna{v)ldV = —Na/V^ = —Ua/V, one recovers 
the compressibility sum rule 



Px ^ = E*^" ~ / '^iT'an^c'^^^{r)dr. 



(3.18) 



Although equation (3.8) is a natural generalization of the compressibility sum rule and similar expressions exist in 
the literature (see for instance Ref. [14]) we included here the derivation of this result in the context of the Poisson- 
Boltzmann and other generic local density functional theories because these theories are approximate (non-exact) and 
nothing garantees in advance the validity of equation (3.8) for non-exact theories. 
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IV. CONCLUSION 

Within the cell model and a generic local density functional theory, we have considered a suspension of charged 

spherical macro-molcciilcs of arbitrary dimension. Wc have cast the corresponding compressibility in a form where the 
sign of this quantity is positive under a (weak) sufficient condition: the grand potential uj appearing in Eq. (2.1) should 
be a convex-up function on densities {no}- This proves that the stability requirement of a positive compressibility is 
fulfilled by Poisson-Boltzmann theory as well as several improvements upon this mean- field approach^^^. Our results 
show that such theories yield stable suspensions, and cannot exhibit a critical point associated with a gas-liquid phase 
separation. 

This result should be put in perspective with the recent proofs that within PB theory, the effective interactions 
between two identical colloids confined in a cylinder of infinite length is necessarily positive^^'^^. This proof has been 
extended to the more general family of local density functional theories (2.1) in Refs. [3,17]. It was shown that the 
positive definiteness of the local free energy density (or equivalently of the grand potential density in the semi-grand 
canonical situation) was a sufficient condition for repulsive interactions. The results derived in this article show that 
under the same circiimstanccs exactly, the bulk modulus within a cell model is also a positive quantity. 

We explicitly considered the semi grand-canonical situation where the macro-ions suspension is dialyzed against an 
electrolyte reservoir. Our results may however be extended to other electrostatic situations, such as the canonical one 
where the mean salt content in the suspension is fixed, or to the situation where the macro-ions are held at constant 
potential rather than constant charge. 

APPENDIX A: ON THE CHOICE OF THE GREEN'S FUNCTION 

We consider the Green's function G(r, r') satisfying 

V?G(r,r') = -5,J(r-r') (Al) 

with yet unspecified boundary conditions. From a standard identity (Green second identity), the solution of Poisson's 
equation VpV = —SvP{y) obeys the relation^ ^ 

V(r)= / p(r')G(r,r')rfr'+i- / ^ll G{r,r') dS' - ^ [ V'(r') |^(r, r') (A2) 

where (?„ = d/dn denotes the normal derivative and dTZ is the surface delimiting the region TZ. 
The Dirichlet Gd and Neumann Gn Green's function satisfy the boundary conditions 

GD{r,r') = and dn' GN{r,r') = -1/ R"'^ ioi r' € dU. (A3) 

For Dirichlet boundary conditions, we therefore have 

V(r) = / p(r') Gn{r, r') dr' ~ ^ [ V(r') ^(r, r') dS' (A4) 

Jtz Jdn '^'"^ 

For the isotropic situation considered here and from Gauss law, the last term in equation (A4) may be written as 

-V(i?)^ / ^(r, r') dS' = ^{R). (A5) 
Jdn c'"' 

So far, the choice Gd{^, r') = for r' G dTZ does not ensure that ijj{R) = 0. We chose to fix the zero of the electric 
potential at the boundary of the cell ip{R) = so that we recover Eq. (2.4). 
For Neumann boundary conditions, 

V'(r) = Ws + [ P(r') G^ir, r') dr' + ^ ( Giv(r, r') dS' , (A6) 

Jn '^•^ Jdn o''^ 

where the first term {tp)g = il){R) is the average of the electric potential on the surface of the cell and we chose it to 
vanish. The last term in equation (A6) is 

^dni^iR) [ GN{v,v')dS' = ( p{v')dv' [ GN{r,r")dS". (A7) 

'-'i/ Jdn -n- Jn Jdn 
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We have used Gauss law: —dnip{R) = J p{v') dr' / R^~^ . A proper choice of the Green's function ensures that the 
surface integral on the right-hand-side of (A7) is independent of r, and vanishes (one can shift Gm by an arbitrary 
constant). Explicitly, in the three dimensional case, the choice 



1 1 °° 

o„(.,.') = -- + -+2; 



£=1 



^ {rr'Y 



Pe{cos0) 



(A8) 



with r> = max(r, r') and r< = min(r, r'), the angle between r and r' and Pi the Legendre polynomial of order £, 
makes the term (A7) to vanish. Finally, with this choice for the Green's function and for the reference potential, we 
have in both cases of boundary conditions 



V(r) 



= / p{r') G{r, r') dr' and tp{R) = 0. 



(A9) 



A useful property that follows from these considerations is that for any isotropic charge distribution p{r) (eventually 
non-globally neutral) we have in both cases of boundary conditions 



p{r')G{K,r')dr' = 



(AlO) 



where R E dTZ. This follows directly in the Dirichlet boundary conditions case from G(R, r') = and in the 
Neumann boundary conditions case form the particular choice (A8). Finally, we emphasize that the symmetry 
property G(r, r') = G(r',r) is not necessarily fulfilled by a generic Green's function^^, but may be imposed as a 
separate requirement, and holds for the functions considered here. 
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